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Abstract 

We consider a Josephson junction with an arbitrary transmission coefficient T) between a singlet 
and a triplet superconductor with the latter order parameter characterized by a d- vector of the form 
(kxi) — kyx). Various quantities such as the tunnelling current, spin accumulation, and spin current 
are calculated via the quasiclassical Green's functions. We also present a symmetry argument on 
the existence of these quantities and their dependencies on the phase difference across the junction. 
A physical picture is also given in terms of the Andreev states near the junction. 

PACS numbers: 74.45. +c,74.20.Rp,72.25.-b 
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I. INTRODUCTION 



Recently, there has been much interest in manipulating the spin degree of freedom of 
electrons in condensed matter systems. Phenomena such as spin current, spin Hall effect, 
spin accumulation and magneto-electric effects have received a lot of attention.- These phe- 
nomena have been discussed in a variety of systems, including metals, semiconductors, and 
even insulators. In this paper, we discuss the spin current and spin accumulation near a 
Josephson junction. We shall in particular consider a junction between an s-wave supercon- 
ductor and a pure triplet superconductor with the latter in the state where the "d- vector" 
specifying the spin structure of the Cooper pairs be given hj d = k^y — kyX. 

We are interested in this kxV — kyX state for a number of reasons. This state corresponds 
to the one-dimensional representation in a crystal with tetragonal symmetry,^- and 
hence is one of the simplest example of a triplet state. This state is also believed to be a limit- 
ing case for the order parameter of the non-centrosymmetric superconductor^ CePt^Pds-^^Si. 
There, due to the absence of inversion symmetry in the normal state, the order parameter is 
believed to be a mixture of s-wave and the p-wave order parameter. [The state k^y — kyX 
is the limiting case (perhaps for small Pd concentration 3 — x) where the s-wave admixture 
is small.] Similar mixing of superconducting order parameter of different parity is also ex- 
pected in compounds such as CeRhSis,'* and in superconductivity found at oxide interface.- 
Spin current generated near the surface of this superconducting state with vacuum have 
been discussed recently by two groups,-"^ with and without the mixing of the s-wave order 
parameter due to the absence of inversion symmetry. 

We generalized these considerations to the case where this superconductor is in contact 
with an s-wave superconductor in the form of a Josephson junction with arbitrary trans- 
mission coefficient V, but for simplicity we shall not include any broken inversion symmetry 
effects in the normal state, hence the bulk superconductors are assumed to be pure singlet 
and pure triplet, respectively. Clearly, in the infinitely high barrier limit, our results would 
just be a special case of Ref.-*^. 

For general transmission however, one expects a proximity effect so that near the interface, 
the system acquires properties of a superconductor with mixed singlet and triplet order 
parameters, similar to the case which arises in non-centrosymmetric superconductors,-!^ 
even though our bulk superconductors are each purely singlet and triplet. Effects that are 
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normally not allowed can now appear due to the lowering of symmetries, similar to the 
electro-magneto effects discussed recently for bulk non-centrosymmetric superconductors.- 
There, in particular, a supercurrent can generate a spin polarization in a perpendicular 
direction. Here, we shall investigate how the spin current (and the spin accumulation) 
depends on (and hence can be manipulated by) the phase difference between these two 
superconductors. 

Our investigation is interesting in another point of view. The state kxV — kyX has two 
counter-propagating edge states of opposite spins near a surface (see below) , in direct analogy 
with the quantum spin Hall state often discussed in the current literature.-*^ Our investi- 
gations here then is analogous to considering an interface between an ordinary "insulator" 
(our s-wave superconductor) and a "quantum spin-Hall insulator" ( our kx]) — kyX supercon- 
ductor). Discussions on this and other related triplet superconductors from this point of 
view can also be found in Ref.— li^. 

A recent paper— also studies the spin accumulation near a Josephson junction between 
a pure singlet and pure triplet superconductor. In that paper, only the very special p- 
wave state where d is independent of the momentum direction k was considered. Spin 
accumulation was shown to exist near the junction, with the spin direction along d. The 
authors suggested the detection of this spin accumulation as a method of identifying triplet 
superconductors. However, the constant d vector is a very special case. A general triplet 
superconductor is expected to have k dependent d vectors.- For these more general cases, it 
is then unclear if spin accumulation would exist, and in which direction the net spin lies. We 
would like to provide a general consideration using this {kxy — kyx) state as an illustrative 
example. 

Our paper is organized as follows. We begin with a symmetry argument in SecHIl We 
then present our calculations with the quasiclassical method in SecJIIII The subsections 
provides our results, first for the special cases of perfect and small transmissions, then the 
more general case with arbitrary T). We summarize in Sec JIVI We employ a generalization 
of the "exploding and decaying trick", which we explain in Appendix lAl 
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II. THE JUNCTION GEOMETRY AND SYMMETRY CONSIDERATIONS 



We shall then consider a Josephson junction between an s-wave superconductor and a 
purely triplet superconductor with d = k^y — kyX. For simplicity we shall consider the two 
dimensional case, or equivalently the three-dimensional case but dispersionless in k^. A 
schematic view of the junction is shown in Fig. [H We shall show that symmetry argument 
forbids existence of certain quantities, and in the case where a quantity is allowed, its 
dependence on the phase difference is constrained. We search for symmetry operations under 
which the junction would map back to itself. Caution has to be made to account for possible 
changes of the phase of the order parameters under these operations. These considerations 
are along the same line as those applied earlier by one of usi^'i^ to the Josephson current 
across a junction. 

The s-wave (triplet) superconductor occupies x < (>)0. The order parameter A is 
a 2 X 2 matrix in spin space. We have, for x < 0, A = As{i(Ty) whereas for x > 0, 
A{k) = Api (d{k) ■ (t^ ay, d{k) = kxH — kyX specifies the triplet structure of the pairs. We 
shall for simplicity ignore anisotropy of the magnitude of the superconducting gaps. In this 
case. As and Ap are independent of k. 

First we consider the time-reversal transformation B under which the supercurrent and 
spin accumulation are odd while the spin current is even. The annihilation operators trans- 
form as Gag^B"^ = and Ga^j^B"^ = — Using the fact that A transforms as the 
corresponding anomalous average, simple algebra then shows that A^ ^ A*, Ap A* with 
d unchanged (using that d is real). Hence the phase difference changes sign. It follows that 
the supercurrent Jj{x) = ~Jj{~x)^ spin accumulation S\x) = ~'S'*(— x), and spin current 
Jjix) = Jji^x) for polarization and flow along i and j, respectively. 

Under a reflection in the x-z plane, the order parameter d = k^y — kyX transforms 
according to {kx,ky,kz) — > {k^, —ky^k^) and (x, z) — > (— x,y, — i), respectively. Hence 
both superconductors are invariant and the phase difference x is also unchanged. The only 
non- vanishing currents, spins and spin-currents allowed are thus Jx,z^ S^, J^^^, and J^'^. 
Since the dispersion in z is not considered, and will not be mentioned hereafter. We 
can also consider a reflection in the x-y plane under which {k^, ky, k^) {kx, ky, —k^) and 
(x, 5) (— X,— z). The resulting order parameter A^ As but Ap ApC*'^, hence 
the phase difference X ^ X + ^- We then have 
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Jxix) = Jxix + tt) , 

j'Ax) = -mx+^), (1) 

Jyix) = -J^iX + ^), 
Jyix) = JyiX + T^) ■ 

Other symmetry operations (such as vr rotation about x) just produce relations that can 
be found by combinations of those hsted above. We note in particular that the spin accu- 
mulation lies entirely along the y direction. In the limit of zero transmission, all quantities 
are independent of x- In this case, all spin accumulations must vanish and the only finite 
spin current is J^. These results hold even when more general components of the order 
parameter (e.g. kxky{kxX — kyy) in Ref.-) are included. As we shall see later, only J^, 
and Jy are found to be finite in our calculations. 

III. QUASICLASSICAL GREEN'S FUNCTION 

We now present our calculations and the quasiclassical method. At positions other than 
the interface, the quasiclassical Green's function g, a function of momentum direction k, 
Matsubara frequency e„ and position r, obeys 

[ienTs-A,g]+tvjCk)-V9 = 0, (2) 
with the normalization condition 



f = -vr^ . (3) 

Here Vf{k) is the Fermi velocity. The boundary condition at a; = will be stated below. 

, A 

A(fc) specifies the off-diagonal pairing field. A = where A isthe 2x2 order 

-A} 

1 \ (00, 

parameter matrix in spin space. With r+ = and r_ = in particle-hole 

oy \^ 1 

space, we then have A = As(zcry)r_,_ + A*(zcr.y)r_ for x < and A = Api{d{k) ■ (T)o-yr^ + 
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A*iay{d{k) ■ (j)r_ for x > 0. In order to have tractable analytic solutions for g, we shall 
also ignore the self-consistent gap equation and hence the spatial dependence of As and 
Ap. We shall also assume for simplicity that the Fermi velocities magnitudes of the two 
superconductors are identical and independent of k. 



A. Perfect Transmission 

In this case the boundary condition at x = is simply that g is continuous. g{k, e„, 0) is 
given by (see Appendix Rl). 

g{k, e„, 0) = -i7r{a, 6}~^[a, b] , (4) 

where [, ] and {, } denote commutators and anti-commutators and d and b the appropriate 
exponentially decaying and increasing solutions along the quasiclassical path. 
For kx > we need d = dp, b = bg in eq. (jl]). We find 

, , d{k) ■ (J 

fi((/c,en,0) = C3(fc,e„)r3 C3(A;,e„) _ + (o.rf.) , (5) 



—ayd{k) 



GO 



y 



where {o.d.) denotes off-diagonal terms in particle-hole space that we would not need, 

if x_ |App|Aspsinxcosx + iien(asap + e^)(ap + as) f„s 
C3[k, en) - -vr -— — — — ,0 „„o , ) W 



'ft \- lA MA I + el)smx + ienjcyp + as)cosx 

where x = Xp ~ Xs the phase difference, as = (e^ -|- lA^^)"^''^, ap = (e^ -|- |App)^^^. The 
result for C3 was also given in Ref.-i^. For k^ < 0, we need d = bp, b = dg in eq. (jl]). 
Alternatively, we can also use the symmetry^^ g{—k, — e„,) = T2^*''(A;, e„)r2 where tr denotes 
the transpose. g{k,en) is still of the form in eq. with C3(— A;, — e„) = — C3(A;, e„) and 
c'^i—k, —e„) = —c'^ik, Cn)- Note that we have defined the C3, C3 coefficients with k dependent 
d vector in eq. ([5]), and d{—k) = —d{k). 

The number current density along x can in general be expressed as 



Jx = ^NfVf J ^(cos0)T^Tr [r3^(A;,e„)] , (8) 
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where is the angle of k with respect to x, Nf is the density of states per unit area for 
a single spin species. The symbol Tr represents taking a full trace in both the spin and 
particle- hole spaces. Only the C3 component in eq. (E]) contributes to Jx- The spin density 
in the i direction at x = can be expressed as^i 
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^^TVTr 
27r ^ 



(9) 



Here we define the symbols a* by a^' = a,j., cyy = cfyT^ and (7^ = 0"^. So here only is finite 
and is associated with C3 in eq. (I7l). The spin current densities, with superscript (subscript) 
denoting the spin (flow) direction at x = 



J! 



'-NfVf 



(10) 



4 27r 

Note that the three components of (TV3 are cr^ra, ay and 0^2X3. It follows that all the 
spin currents vanish since g of eq. ([5]) does not contain any aVs components. Physically, 
the Andreev equation for each k is decoupled from other paths, and hence can be block- 
diagonalized using quantization axis along d{k). Along this axis, both the singlet and triplet 
superconductors consist of only ][ pairs. These Cooper pairs do not have any net spins, and 
they cannot contribute to any dissipationless spin current. See also the discussions near the 
end of subsection C. 

Next we present explicit results for the case of equal gaps on both sides, i.e. [A^l = 



I A, 



|A|. Here the interface bound states, which correspond to the poles of g in eq. 



([5]), are essential for the quantities in eq. ([8]) and ([9]). It can be shown that for the right 
moving path [k^ > 0), the bound states of spin parallel and antiparallel with d{k) are given 



by Eb 



I A| cos(|-)sgn 



sm 



and Efj^i = | A| sin(|-)sgn 



cos 



left moving path [k^ < 0), the bound state energies are 



respectively. For the 
and 



I A| sin(|)sgn 



cos(|) 



E 



'b,l 



|A| cos(|)sgn 



sin(f) 



. Notice that we adopt a common spin quantization axis for 
both right and left moving paths (caption of Fig. [2]) to facilitate the following discussions. 
The bound state spectra are plotted as a function of phase difference x Fig- El It can be 
seen that, for a given path, the two branches of opposite spin projections are identical except 
separated by tt, which reflects the invariance of triplet order parameter under Xp ^ Xp + ^ 
and d —d. 

The analytical results for is obtained using eq. ([8]), which gives 
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2|A| vr 



cos(x/2) tanh _ sin(x/2) tanh ^^^^^Z^)' 



'111 



2T ' 2T 

where i?Ar denotes the corresponding resistance in the normal state. Eq. (fTTj) coincides with 

the previous results in Ref.-^^. Jx is plotted in Fig. [3] and the present case corresponds to 
the line denoted hj V = 1. can be understood by summing over contributions |^ from 
occupied bound states. Notice that a current jump occurs whenever x is a multiple of vr. 
When X is slightly larger than 0, the state labelled the red square in the left panel and the 
one labelled by the green circle in the right are occupied. Only the latter bound state with 
a negative slope contributes to J^.. When x is slightly less than 0, on the other hand, the 
black square in the left panel with a positive slope is occupied and contributes to J^- 

Moreover, the splitting between bound states actually contributes to a finite spin accu- 
mulation near the interface along some direction. Consider < x < and zero temperature. 
Referring to Fig. [21 for the right and left moving paths, the states with spin parallel to the 
quantization axis defined in the caption are both populated. As the parameter varies be- 
tween ±7r/2, this quantization axis varies. A net spin is generated along the positive y-axis, 
whereas the x component adds to zero. Analytically, the spin accumulation can be obtained 
from eq. 0, which gives 



sy = hNf\A 



cos(x/2) tanh ^^M^ + sin(x/2) tanh ^^^^^Z^) 



;i2) 



2T ' 2T 

As a function of x, the spin accumulation S"^ for both sides of the interface is plotted in Fig. HI 

The present case corresponds to the line with T> = 1. In addition, 5*^ is also continuous across 
the interface for perfect transmission. We note however that, if the magnitude of the gaps 
of the two superconductors are unequal, there can also be contributions due to continuum 
states, as in the case of supercurrent between two unequal gap s-wave superconductors.— 
Since the Green's function decays as e-2akl/-«/|cos0|^ (jgc^ys in a distance of order of 
coherence length hvf/\A\ away from the interface. The total spin accumulation is of order 
h^NfVf per unit length along the junction. 



B. No transmission 



In this case all particles are reflected. The behavior of the s-wave superconductor for 
X < is trivial and we shall thus concentrate only on the triplet superconductor on the 
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right. Let us denote the incoming wavevectors by k and the reflected outgoing wavevectors 
by k, with A;^; > and k^ < 0. See Fig. [H We label positions along the quasiparticle path 
consisting of each pairs of k and k by u, with u < {u > 0) labels the part for k (k). g{u) is 
continuous at m = 0, and can be obtained from eq. (jlj) with a <x.p{k) and h bp(k). Since 
d{k) 7^ d{k), we shall introduce the quantities C = d{k) ■ d{k_) and D = d{k) x d{k). Note 
that + |Dp = 1. The part of ^(0) which is diagonal in particle-hole space and even in 
is found to be 

|ApP ((D-a) \ 

''[2el + \Apmi + C)]y ay{D.a)ay)' 

For our stcitB, C — — cos20 and D = zsin20. It follows that there are no currents Jj. For 
a given pair of wavevectors k and k, there is in general a finite spin along D \\ z. However, 
the contribution from the pairs of wavevectors in opposite directions sum to zero (That is, 
between the pair with outgoing k and incoming wavevector being —k, or alternatively, ±0). 
Therefore 5*^ = 0, and there is no spin accumulation in any direction, which can also be 
seen by noting that g does not contain any a* component. The only finite spin current is Jy 
associated with the a^r^ component, and its value at a; = is given by 

- T '-'■'^'^?- |2.;+'|A'ni + c)] • 

where is the z component of D. Since g(k, e„, 0) = g(k, e„, 0), the angular integral in eq 
(fT3|) has been replaced by twice the contribution due to outgoing wavevectors. The factor 
sin0 is due to ky = ky At zero temperature, the spin current density Jy = fiNfVf^A^ at 
the interface and decays into the bulk within a coherence length. The total spin current is 
of order h^Njv'j. 

The physical picture of the spin-current is similar to that of the edge current in the so- 
called chiral superconductors,— and has been discussed also in, e.g., Ref.-. Our triplet state 
consists of It pairs and || pairs only, with wavefunctions respectively given by (— rf^'+irf^) — »• 
i{kx — iky) and {d^ + id^) i{kx + iky). Due to the phase difference between the order 
parameters of the incoming and outgoing momenta, each spin component has a bound 
state (for a given pair of incident and reflected wave-wavevectors) but opposite energies: 
e = =F|Ap| sin0 for spin up (down) respectively near the surface. Thus the up (down) spins 
preferentially occupy the states with positive (negative) y momentum, contributing to a net 
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spin current Jy. In this picture, it also follows that Jy vanishes for V = 0. 
C. General Transmission 

In this subsection we consider a general interface between our singlet and triplet super- 
conductor of (angular and spin independent) transmission coefficient V. We denote the 
incoming (outgoing) wavevector on the right by k and k, and conversely for the left, see 
Fig. [TJ The corresponding Green's functions g{k, x = 0±), g{k, x = 0±) on the two sides of 
the spin-inactive interface are related to each other by a set of boundary conditions given in 
Ref.— . See also Appendix |Xl It is more convenient to express them in terms of the difference 
gd = g{k,x = 0+) — g{k,x = 0+) = g{k,x = 0_) — g{k,x = 0_), which is continuous across 
the interface, and the sums s"^^^^ = g{k,x = 0_|_(_)) + g{k,x = 0+(_)). It can be shown that 
the supercurrent Jx in eq. ([H]) and spin currents J* in eq. (|TOl) across the interface can be 
expressed solely in terms of the difference cjd- Note that which specifies the angle for k 
is now restricted within ±7r/2. Moreover, the and r^a'' components of cjd are associated 
with Jx and J*, respectively. Below we consider the equal gap case for simplicity in which 
c/d can be worked out analytically via eq. flA16p . The rs component contributing to Jx is 
given by 





( 


-7r)p2 


A 


^sin(2x) 




4sin2x + 4(l -r')a2 


A 


^ sin^ 



By numerically performing the sum over the Matsubara frequencies, Jx for arbitrary V is 
plotted in Fig. [31 Note that the current is odd and is periodic in the phase difference x with 
period vr, as noted also in Ref.—. See also SecHIl Second, we find that none of the r^a^ 
components appear in gd, and hence all the spin currents J* across the junction are zero. 
We note, however, that the (spatial) symmetry argument in Sec. [ITl allows a nonzero J^ as in 
eq. ([1]). Therefore, the vanishing of J| results from other symmetries which we shall discuss 
near the end of this section. 

At the right side of the interface, the spin accumulation 5* and the spin current J* fiowing 
parallel to the interface can all be expressed in terms of solely. Here the a* components 
in s*" are needed for and the r-^cr^ ones are for J* as required in eq. Qj and ffTOj) . By using 
eq. (lA17p . the a* components are listed below. 
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—i sin (j) cos x<^^n 




+ 


a - f + f 6^; 


cos sin X o-yT^ 




A 


4sin2x + 4(l 


\A 


^ sin^ (p 



The spin accumulation 5*^ is identically zero because the coefficient in ax is odd in e„, and 
the factor factor sin also gives zero after the angular integration. This result is consistent 
with our symmetry argument in Sec [Tll The only finite spin accumulation is which is 
shown in Fig. H] due to the ayT^ component in eq. (fTSll . Note that S'^(x) obeys the symmetry 
in Sec [Til and has period 2tt. As for the spin current, the only nonvanishing component of 
TaO"* is given by, 



For V = 0, Jy does not depend on x- For I' < 1, the phase dependence comes from the 
sin^ X term in the denominator. The Jy versus the phase difference x is plotted in Fig. [5] for 
various V. This spin current is even in x and is periodic with period vr. (see Sec |Tl])- We 
note that the vanishing of r^a^ components leads to zero Jy which was not anticipated by 
our symmetry argument in Sec. HIl 

At the left side of interface, 5** and J* can be calculated via in eq. ( 1A18I) . The a* 
components are identical to those in except that a and e„ are interchanged in the square 
bracket [...] of eq. (fT5|) associated with CyT^ component. The numerical results for S^{x = 0_) 
are also shown in Fig. HI It can be seen that is continuous across the interface only for 
V = 1. In addition, all the terms associated with the spin current Jy, including the cTzTs 
component, vanish for all V. Consequently, all the spin currents vanish on the left side. 

The vanishing of Jj for a; < and J* for all a; is a result of spin-conservation. Observing 
that (TajTa, ay, a^T-i commute with rs and A^, we see that, by multiplying eq. ([2]) by these 
matrices and then taking the trace, Vf ■ V(Tr[(jV3^]) = 0. That is, the spin current is 
constant along any quasclassical path at any point inside the singlet superconductor. Since 
the spin-current vanishes at x — * —00, it follows that the spin current on each quasiclassical 
path vanishes for x < 0. Hence Jj = for alH, j if x < 0. Note that this vanishing of 
the spin-current does not rely on angular integration. Since J* = for x = 0_ and the 
spin current is continuous across an spin-inactive interface {gd is continuous), J* = also 
for X = 0+. At any point x > 0, the Green's function is a linear combination of its value at 
X = 0+ and x ^ 00 where J* also vanishes. Hence = also for all x > 0. The vanishing 
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of Jj at a; = 0_ and J* at a; = can also be easily proven using eq. (1A16P and (lAlSP using 
that fact that aVs commutes with gl^^ and TilaW^gl^^] = 0. 

As mentioned, the symmetry allowed Jy is found to vanish in our calculation. We have 
checked that the vanishing of is also true in the case of \As\ ^ |Ap|. We do yet not have 
a simple physical explanation of this result. Mathematically, this follows from the fact that 
absence of the cTxTs term for g on the left of the interface (due to spin conservation) is carried 
over to g on the right. Vorontsov et alA have also considered the interface between vaccum 
and a noncentrosymmetric superconductor with finite spin-orbital Rashba energy, which 
lifts the energy degeneracy between quasiparticles at the same momentum but opposite spin 
projections. They showed that this can lead to some finite and oscillating Jy and J^. We 
expect that this may also happen in our junction. 



IV. DISCUSSIONS AND CONCLUSIONS 



We have considered the spin accumulation and spin current near a Josephson junction 
between a singlet and triplet superconductor. We showed that symmetry arguments (Sec. [IT|) 
place strong restrictions on the existence of above physical quantities and their dependence 
on phase difference x across the Josephson junction. Comparing with the pervious work,— 
this method also applies for any triplet pairing wavefunction and provides a more general 
way of determining the direction in which the spin lies. Conversely, the direction and phase 
dependence of the spin accumulation actually inform us about which symmetry is broken 
by the junction and hence the symmetry of the triplet order parameter itself. Moreover, 
the quasiclassical Green's function technique is employed to quantitatively investigate the 
predicted supercurrent J^, spin accumulation 5*^, and spin current Jy'^. Jy turns out to be 
zero for our junction, though it is symmetry allowed. For transmission coefficient < < 1 
in our calculation, the spin accumulation 5*^ and spin current Jy coexist within a coherence 
length at the triplet side, a feature which does not appear in the previous studies.-*^ 

In conclusions, we have calculated the spin accumulation and spin current near the inter- 
face of a singlet-triplet junction with the triplet order parameter specified hj d = k^y — kyX. 
The method of quasiclassical Green's functions as well as the symmetry arguments can be 
generalized to other junction with arbitrary pairing symmetries. These spin accumulation 
and dissipationless spin currents depends on the phase difference and hence can be controlled 
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by the charge current passing through the junction. 



V. ACKNOWLEDGEMENT 

This research was supported by the National Science Council of Taiwan under grant 
number NSC-95-2112-M001-054-MY3. 

APPENDIX A: EXPLODING & DECAYING TRICK 

In this Appendix we explain the exploding and decaying trick. This trick has been used for 
pure s-wavei^ and pure p-wave pairing (e.g., Ref.— ). From these references, one can actually 
deduced that the method can be generalized to mixed singlet and triplet pairs, so that results 
such as eq. can still be used. However, we would like to provide our alternate derivation 
below to show that it is indeed applicable for mixed pairing, and moreover we believe that 
our presentation may be more transparent to some readers than those in the literature. 
We also note that this method is not limited to spatial independent gaps, though we shall 
discuss only the (piecewise) constant gaps case to simplify our presentation. Furthermore, 
this method can be easily implemented numerically, as has been performed in, e.g.,—^^^^ 
etc. 

We begin by reviewing the first the trick for pure s-wave superconductor. Writing u as 
the parameter along a quasiclassical path, eq. ([2]) can be written as 

[ienTs- As{iay)r+- A*{iay)T_, g{u)]+ivfdugiu) = (Al) 

where we have suppressed the k and e„ dependence of g. A "constant" solution (satisfying 
also eq. ([3])), which is also the g for a bulk superconductor, is given by 



_2e„r3 - As{iay)T+ - A*{iay)T^ 
(4 + |A. 



9sMlk - , n . , . 12x1/2 y^^) 



and is thus a linear combination of r^jCXyT^ and ffyT^ matrices only. It is also possible to 
find other solutions to eq. ([2]) (without satisfying eq. ([3])) which are linear combination of 
these three matrices only. They are, with = (e^ + lA^p)^''^, 
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a,{u) = e-2"»"/"^' (-z|A,|V3 - A,(a, + e„)«(Tj,r+ + A^a^ - e„)2(T^r_) (A3) 
bsiu) = e+2-=«/-/ (+z|A,|V3 - A,(a, - e„)z(T,r+ + A:(a, + e„>a,r_) (A4) 

which will be called the decaying and exploding solutions "in the same block"—. We note 
that they satisfy =, P = 0, {^s,buik) = {fi's.buik, b} = 0. In fact, ^s,buik can be written as 
g = —m{Pi — P'l) with P\ = db/{d,b} and P2 = ba/{d,b} being projection operators with 
P^ + P^ = 1, P1P2 = P2P1 = 0, and aA = 0, dp2 = d, bPi = b, 6P2 = (see, e.g., Ref.-^). 
Similar results apply to the pure triplet superconductor. The bulk solution is 

— * — * 

ienT3 - Ap{id ■ aay)T+ - A*Jiayd ■ a)T_ 
c/p,buik = -vr , , . , , ,1/2 (A5) 



{el + I A 



pi 



and is thus a linear combination of rs, {d ■ a)ayTj^ and (yy{d ■ a)T^ matrices only. The other 
solutions to eq. ([2]) (without satisfying eq. ([3])) which are linear combination of these same 



three matrices are, with Up = (e^ + |Ap 



2\l/2. 



dp{u) = e 2"''"/^^ (-i|Ap|V3 - Ap{ap + ■ o)oyT^ + A;(ap - er,)iay{d- a)r_)(A6) 

bp{u) = e+2"p«/-/ (^+^|Ap|V3 - Ap{ap - e^W ■ a)ayT+ + A^ap + en)iay{d ■ a)r_)(A7) 

Let us now consider our junction, and begin with the case of perfect transmission. A = 
As^p for X < (>)0, and g is continuous at x = 0. Let us first consider kx > 0, and label a 
point on the quasiclassical path by u, with m = at the interface, (hence u = x/kx). g must 
decay to ^s,buik (s'p.buik) as m — —00 (+0x3). We note however, that we cannot just try the 
ansatz g{u) = gp^huik + Cpdp{u) for -u > and g{u) = g s, hulk + Cbds{u) for m < for some scalar 
coefficients Cg and Cf,. This is because the matrices involved for u > (<)0 are then different, 
so g being continuous at u = can never be satisfied. To explain more clearly our idea of 
solving this problem, let us first consider the special case d = z, so that i{d ■ a)ay = a^. 
Then 5'p,buik, dp above are linear combinations of t^jCXxT^. To find a possible continuous g 
at M = 0, we must therefore include also decaying solutions for u > which also involves 
(JyT± (due to the singlet superconductor on a; < 0), and exponentially increasing solution 
for M < which also involves axT±. One can find these solutions easily, as done explicitly 
in ReL-^. We can however also note that these needed solutions can be written as (cT^ra)^^ 
and {(JzT3)bs- (Note that (c^ra) commutes with ts, A^ and Ap). Hence we can tiy^ 
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g{u) = ^s.buik + Csbs + (s{(TzT3)bs u<0 (A8) 
g{u) = gp,hu\k + Cpdp + Cp{cr^T3)dp u> (A9) 

where Cg^p and (s,p are scalar coefficients to be determined. Note that now g for both u < (>)0 
consist of Ts, CT3, <JyT± and cr^r± matrices and hence a solution is possible. Note that eq. ([3]) 
is satisfied. Since ^(0) can be expressed as either eq. ( 1A8I) or ( ]A9I) . we can determine the 
coefficients c^^p and (s,p using simple algebra, but a simpler procedure is to left-multiply eq. 
(lASP and (]A9p (at m = 0) by 63 and respectively to obtain 



bsg{0) = 6s^s,buik = -inbs (AlO) 
apfi'(O) = ap^p,buik = +«7rap . (AH) 

Note that the unknown scalar coefficients have all disappeared. Further multiplying eq. 
(lAlOp and (lAlip respectively by and bp, and adding the two equations, we obtain thus 



g{0) = -^^{apAy'iapX] ■ (A12) 

Repeating the above procedure by post- rather than pre- multiplication actually shows that 
we can also reverse the order of the anticommutator and commutators in eq. (]A12|) . as can 
be also verified explicitly. Note that {dp,&s} is a linear combination of 1 and a^r-^ only. 

For kx < 0, u < 0(> 0) corresponds to x > (< 0). Following again the above procedure 
and ensuring that the solutions decay correctly to their respective bulk values at n ^ =pcxo 
gives us the analogous formula 

^(0) = -iTr{ds,bpy^[ds,bp] . (A13) 

Eq. (]A12I) and (1A13I) are the special examples of eq. (jlj) in the present case. For general d{k), 
to ensure the continuity of 5^ at x = 0, we need matrices T3, cryT±, {d-a)o-yT^ and ay{d-a)T^. 

. . I (d-a) \ 

A matrix that commutes with t^, Ag, Ap can be seen to be ^ = Si. 

V ay{d-a)(yyj 

The argument above can be repeated with this matrix replacing cTzTs above. 

The above argument actually does not depend on the fact that cr^rs (or Ei defined 
above) be common to both sides of eq. (lASp or (1A9p . To see this, let us ffist consider the 
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singlet superconductor. We note that the matrices 1, ay, a^T^,, cTzT-s all commute with t^, 
CTyTj-, SO they are automatically solutions to eq. flAll) . Since the product of two solutions 
to eq. (jAl]) is also a solution, we see that ^s,buik, o-j/5's,buik, (^x^sgsMik, '^zisS's.buik are also 
"constant solutions". There are also in fact four decaying solutions dg, aydg, axTsds, cTzT^dg 
and similarly four exponentially increasing solutions. Note that we now have 16 solutions to 
the 4x4 matrix equation flAip . and hence any solution to eq. (lAip can be written in terms 
of them. The most general g which decays to ^s,buik at u ^ —oo (for kx > 0) can be seen to 
be 

g{u) = gs,hulk + Csk + CsACTyk + Cs,20-xT3bs + CsO^^Tsfe^ (A14) 

where Cg, Ci-3 are scalar coefficients. Note that no constant solution other than ^s,buik can 

appear on the right hand side of eq. flA14p due to the condition at u — oo. Since I, (Jy, 

(JxT^, a^Tz all commute with rs, ayT±, they commute with hg. Left multiplication of eq. (1A14P 

with hg again yields eq. ( lAlOp . 

The triplet superconductor on x > can be treated similarly. For a given A;, we have 

^ [d(k) ■ a) 



already noted that the matrix 



Si (A;) commutes with rs, 



ay{d{k) ■ a)ay 

{d{k) ■ a)ayT^ and ay{d{k) ■ a)T_. Two other matrices with this property are (besides 1) 

! (^2,3 ■ ff) \ . . 

= ll2,'i{k) where ^2,3 are the two vectors orthogonal to d{k). 

^ -cry{d2,3 ■ ^)(^y J 

These four matrices 1, Si 2,3 (A;) are trivial solutions to eq. ([2]) for the triplet superconductor. 
Four other constant solutions are the product between them and ^p,buik- Again there are 
four decaying (increasing) solutions obtained by their product with dp (bp). We again have 
a total of 16 solutions to eq. ([2]) for the triplet superconductor. The most general solution 
to g{u) with g{u) ^p,buik as m — > 00 (again for kx > 0) is 



g{u) = ^p,buik + Cpdp + Cp,i^i{k)dp + Cp^2'^2{k)dp + Cp_3S3(A;)ap . (A15) 

On noting that Si 2.3(A;) commute with r^, {d{k) ■ a)ayT^ and ay{d{k) ■ (t)t^ and hence dp, 
left multiplying eq. (lAlSp by dp again yields eq. flAlip . The rest of the demonstration of eq. 
(1A12P goes through unchanged. Similar argument applies for kx < 0. 

For finite transmission V, g is now in general discontinuous at x = 0, and g for incoming, 
refiected and transmitted paths are all related. A general boundary condition non-linear 
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in g was first derived independently by Zaitsev^ and Kieselmann^^. In Ref.— a simplified 
linearized form of the boundary condition was provided. The derivation given there was 
for singlet superconductors. However, the boundary conditions derived by Ref.— were 
actually independent of the assumption on the parities of the superconductors. This can also 
be checked by using the formulas derived by Millis et al?'^ for a spin-active interface between 
two superconductors of different parities. By ignoring the spin dependence of the scattering 
amplitudes in Ref.— and eliminating the "drone amplitudes" there, one can show that the 
nonlinear boundary condition of Kieselmann^^ can be recovered. This non-linear boundary 
condition can then be linearized using arguments similar to those used in Ref.—: We express 
g in the form eq. flA14p and (lAlSp for each of the quasiclassical incident, reflected and 
transmitted path but with ^buik replaced by ^aux, the "auxiliary" solution corresponding to 
the completely reflecting case^° (That is, for example, ^^ux solves the quasiclassical equation 
on the quasiclassical path formed by k_ and k with the physical order parameter on the 
right but with the boundary condition ^aux(^) = .^Lx(Al) at x = 0+). The decaying and 
exploding terms can be eliminated using projection operators^^ with arguments similar to 
those explained above for the perfect transmission case. Thus the derivation in Ref.= can 
be carried over to our present situation. It is most convenient to write the final results in 
terms of s*"'' = g{k, 0±) +g{h-, 0±) at x = 0± and gd = g{k, 0±) —g{k, 0±) where k (k) denotes 
outgoing reflected (incoming incident) wavevector on the right (r). k (k) is also the incoming 
(reflected) wavevector on the left (/). gd (denoted by d in Ref.—) is continuous across the 
interface and is given by 



471-2 vi/aux Jaux/ 



whereas 

[2 - i^)^Lx + ^gi 

and 



_ (2-p)^Lx+p^ :ux 

4jf2 Vi/aux i/aux/ 



, ^g 

aux /■ A 1 o N 



where the subscripts "aux" denote the solution to the P = problem. Since ^aux commute 
with the anti-commutator {^aux5 .^aux}; need not specify the relative order between the 
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numerator and the denominator in eqs. (1A16P - (1A18|) . 

By some straightforward algebra, the complete quasiclassical Green's function gl^^ for 
V = problem in Sec. IIIIBl is shown to be, 



el + |AJ2sin^ 



-[apenTs + -|Ap|^sin(20)cr^r3 



(A19) 



+ Ap(ap sin Ca; + e„ cos (p (7y)ayT^ + A*cry (ctp sin (p ax — en cos crj,)r_ 



which can be shown to satisfy (^aux)^ = — vr^. Together with the trivial ^^ux = -^{.^v.T'i — 
/^s^yT^ — A*ayT^) for the left side, one can obtain ga in the following form. 



g, = {-tnV)CA 



where the matrix C is the inverse of 



(7-1 = (2 - [e2 + (|Ap|sin 



V 



-B 



The matrix B comes from the anticommutator and is given by. 



B 



4 + (l^p|sin> 



{^aux' ^aux} 



-TT" 



= 2Q;pe^ + ie„| App sin(20) cr^ — 2ie,„ | A <, 1 1 Ap | cos sin x 

- 2ap|A5||Ap| sin cos X a^Tz + | App sin(20)(A* cr^r_ - A^ a^r^ 



A is from the following commutator. 



(A20) 



(A21) 



(A22) 



A = a. 



4 + d^plsin' 



S'aux' fi'aux 



= — 2e„|As||Ap| cos cos X CTyVs — 2ia;p| As| | Ap| sin sin x 
- 2 [e„Ap(ap sin + e„ cos ay) + ape„A^] ayr+ 



— 2a,, 



e„Ap(-ap sin + e„ cos ay) - a;pe„A' 



(A23) 



Similarly, the sum can be expressed as. 



4C 



(-^ 2 )^aux 2 9a\ix 



"s(e^ + I Ay 



sm 



(A24) 
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and the expression for s'' is identical to the above with interchange of (1 — ^) and ^ in the 
bracket. 
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FIG. 1: A schematic view of the singlet-triplet junction. The triplet superconductor, with an 
order parameter of magnitude |Ap| and phase Xpi occupies the right (x>0) while the singlet one, 
whose respective values denoted by |As| and Xs, occupies the left (x<0). The quasiclassical path 
is denoted by the direction of quasiparticle momentum k. The angle is defined with respect to 
the X-axis. Incoming and outgoing paths labelled by k and k, respectively, are used for interface 
with non-perfect transmission. 
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FIG. 2: (Color online) The interface bound states associated with the kx > and kx < paths in 
the perfect transmission case. To facilitate the discussions for spin accumulation and spin current, 
I spin means parallel to d{k) associated with the right moving path, i.e. kx > 0, but antiparallel 
to d{k) if kx < 0. 
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FIG. 3: (Color online) The calculated supercurrent Jx for various transmission coefficients P. The 
gap on both sides are set to equal to |A| and T = [A|/100. 
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FIG. 4: (Color online) The spin accumulation at x = O^z on the two sides of the interface for 
various transmission coefficients T> with T = |A|/100. 
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FIG. 5: (Color online) The spin current at x = 0+. T=|A|/100 here. On the singlet side of 
interface, the spin current is zero for all T>. 
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